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Abstract 

A MIMO radar system is proposed for obtaining angle and Doppler information on potential targets. Transmitters 
and receivers are nodes of a small scale wireless network and are assumed to be randomly scattered on a disk. The 
transmit nodes transmit uncorrelated waveforms. Each receive node applies compressive sampling to the received 
signal to obtain a small number of samples, which the node subsequently forwards to a fusion center. Assuming 
that the targets are sparsely located in the angle-Doppler space, based on the samples forwarded by the receive 
nodes the fusion center formulates an £i -optimization problem, the solution of which yields target angle and Doppler 
information. The proposed approach achieves the superior resolution of MIMO radar with far fewer samples than 
required by other approaches. This implies power savings during the communication phase between the receive nodes 
and the fusion center. Performance in the presence of a jammer is analyzed for the case of slowly moving targets. 
Issues related to forming the basis matrix that spans the angle-Doppler space, and for selecting a grid for that space 
are discussed. Extensive simulation results are provided to demonstrate the performance of the proposed approach at 
difference jammer and noise levels. 

Keywords: Compressive sampling, MIMO Radar, DOA estimation, Doppler estimation 

I. Introduction 

Multiple-input multiple-output (MIMO) radar systems have received considerable recent attention, e.g., [l]-[3]. 
Unlike a conventional transmit beamforming radar system that uses highly correlated waveforms, a MIMO radar 
system transmits multiple independent waveforms via its antennas. A MIMO radar system is advantageous in two 
different scenarios [4] -[6]. In the first one [4], the transmit antennas are located far apart from each other relative to 
their distance to the target. This enables the radar to view the target from different directions simultaneously. The 
MIMO radar system transmits independent probing signals from decorrelated transmitters through different paths, 
and thus each target return carries independent information about the target. Combining these independent target 
returns results in a diversity gain, which enables the MIMO radar system to reduce target radar cross section (RCS) 
scintillations and achieve high target resolution. In the second scenario [5], a MIMO radar is equipped with M t 
transmit and N r receive antennas that are close to each other relative to the target, so that the RCS does not vary 

'Copyright ©2008 IEEE. Personal use of this material is permitted. However, permission to use this material for any other purposes must 
be obtained from the IEEE by sending a request to pubs-permissions@ieee.org. 

This work was supported by the by the Office of Naval Research under Grant ONR-N-00014-07-1-0500, and the National Science Foundation 
under Grants CNS-06-25637 and CNS-04-35052 
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between the different paths. In this scenario, the phase differences induced by transmit and receive antennas can 
be exploited to form a long virtual array with M t N r elements. This enables the MIMO radar system to achieve 
superior spatial resolution as compared to a traditional radar system. MIMO radar can achieve a desired beampattern 
by transmitting correlated waveforms [7] -[9]. This is useful in cases where the radar system wishes to avoid certain 
directions because they either correspond to eavesdroppers, or are known to be of no interest. In this paper we 
consider closely spaced transmit and receive antennas and uncorrected transmit waveforms. 

Compressive sampling (CS) [10]-[12] has received considerable attention recently, and has been applied suc- 
cessfully in diverse fields, e.g., image processing [14] and wireless communications [15][16]. The theory of CS 
states that a X-sparse signal x of length TV can be recovered exactly with high probability from 0(K log N) 
measurements via l\ -optimization. Let $ denote the basis matrix that spans this sparse space, and let $ denote a 
measurement matrix. The convex optimization problem arising from CS is formulated as follows 

min||s||i, s.t. to y = $x = $\frs (1) 

where s is a sparse vector with K principal elements and the remaining elements can be ignored; <I> is an M x N 
matrix with M <C N, that is incoherent with 'J. 

The application of compressive sampling to a radar system was recently investigated in [17], [18] and [19]. 
In [17], in the context of radar imaging, compressive sampling was shown to have the potential to reduce the 
typically required sampling rate and even render matched filtering unnecessary. In [18], a CS-based data acquisition 
and imaging algorithm for ground penetrating radar was proposed to exploit the sparsity of targets in the spatial 
dimension. The approach of [18] was shown to require fewer measurements than standard backprojection methods. 
In [19], CS was applied in a radar system with a small number of targets, exploiting target sparseness in the time- 
frequency shift plane. The work of [20] considered direction of arrival (DOA) estimation of signal sources using 
CS. Although [20] focussed on communication systems, the proposed approach can be straightforwardly extended 
to radar systems. In [20], the basis matrix if? was formed by the discretization of the angle space. The source 
signals were assumed to be unknown, and an approximate version of the basis matrix was obtained based on the 
signal received by a reference vector. The signal at the reference sensor would have to be sampled at a very high 
rate in order to construct a good basis matrix. 

In this paper, we consider a small scale network that acts as a MIMO radar system. Each node is equipped with 
one antenna, and the nodes are distributed at random on a disk of a certain radius. Without any fixed infrastructure, 
the distributed antennas in this small network render such MIMO radar more flexible than a fixed antenna array 
since we can choose the nodes freely. For example, the network nodes could be soldiers that carry antennas on 
their backpacks. We refer to such a MIMO radar system a distributed MIMO radar. The nodes transmit independent 
waveforms. We extend the idea of [20] to the problem of angle-Doppler estimation for MIMO radar. Since the 
number of targets is typically smaller than the number of snapshots that can be obtained, angle-Doppler estimation 
can be formulated as that of recovery of a sparse vector using CS. Unlike the scenario considered in [20], in MIMO 
radar the transmitted waveforms are known at each receive node. This information, and also information on the 



November 25, 2009 



DRAFT 



4 



location of transmit nodes, if available, enables each receive node to construct the basis matrix locally, without 
knowledge of the received signal at a reference sensor or any other antenna. In cases in which the receive nodes 
do not have location information about the transmitters, or they do not have the computational power, or they 
face significant interference, the received samples are transmitted to a fusion center which has access to location 
information and also to computational power. Based on the received data, the fusion center formulates an augmented 
^-optimization problem the solution of which provides target angle and Doppler information. The performance of 
l\ -optimization depends on the noise level. A potential jammer would act as noise, and thus affect performance. We 
provide analytical expressions for the average signal-to-jammer ratio (SJR) and propose a modified measurement 
matrix that improves the SJR. For the case of stationary targets, the proposed approach is compared to existing 
methods, such as the Capon, amplitude and phase estimation (APES), generalized likelihood ratio test (GLRT) [2] 
and multiple signal classification (MUSIC) methods, while for moving targets, comparison to the matched filter 
method [21] is conducted. 

Preliminary results of our work were published in [22]. Independently derived results for MIMO radar using 
compressive sampling were also published in the same proceedings [23]. The difference between our work and [23] 
is that in [23] a uniform linear array was considered as a transmit and receive antenna configuration, while in our 
work we focus on randomly placed transmit and receive antennas, i.e., an infrastructure-less MIMO radar system. 
Further we study the effects of a jammer on estimation performance. 

The paper is organized as follows. In Section II we provide the signal model of a distributed MIMO radar system. 
In Section III, the proposed approach for angle-Doppler estimation is presented. In Section IV we derive the average 
SJR for the proposed approach and also discuss a modification of the random measurement matrix that can further 
improve the SJR. Simulation results are given in Section V for the cases of stationary targets and moving targets. 
Finally, we make some concluding remarks in Section VI. 

Notation: Lower case and capital letters in bold denote respectively vectors and matrices. The expectation of a 
random variable is denoted by E{-}. The superscript (•) and Tr( ) denote respectively the Hermitian transpose 
and trace of a matrix. 

II. Signal Model for MIMO Radar 

We consider a MIMO radar system with M t transmit nodes and N r receive nodes that are uniformly distributed 
on a disk of a small radius r. This particular assumption will be used in Section IV for the analytical evaluation 
of the proposed approach. For simplicity, we assume that targets and nodes lie on the same plane and we consider 
a clutter-free environment. Perfect synchronization and localization of nodes is also assumed. The extension to the 
case in which targets and nodes lie in 3-dimension space is straightforward. Let (r|,a*) and (r[,a[) denote the 
locations in polar coordinates of the i-th transmit and receive antenna, respectively. Then the probability density 
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t It t It 

functions of r-' and a- are 

o t/r 

£ t t/r\ ^'i f\ ^ t/r „ 

f r t/r{rJ ) = 0<r/ <r 

and f a t/r(al /r ) = — , -tt < a* /r < tt. (2) 

Let us assume that there are K point targets present. The fc-th target is at azimuth angle 9 k and moves with 
constant radial speed v k . Its range equals dk{t) = (4(0) — v k t, where d k (0) is the distance between this target and 
the origin at time equal to zero. Under the far-field assumption, i.e., dk(t) ^> ?**^ r , the distance between the i\h 
transmit/receive antenna and the fc-th target d\ k ld r ik can be approximated as 

4l r (t) w d fc (t) - vt /r (0 k ) = 4(0) - v k t - 4 /r (e k ) (3) 

where 7] t / r (6 k ) = r\' T cos(0 k - aj r ). 

Let Xi(i)e j2,r ^* denote the continuous-time waveform transmitted by the z-th transmit antenna, where / is the 
carrier frequency; we assume that all transmit nodes use the same carrier frequency and also that the Xi(t) is 
periodic with period T and narrowband. Besides, we also assume the slowly moving targets, i.e., ^-Cl. 

The received signal at the fc-th target equals 

Mt d l (t) 

Vk(t) - /3 fc V^(t-d* fe (t)/c)exp(j27r/(i-^i)), k = l,...,K (4) 

U 

where {(3 k , k = 1, . . . , K} are complex amplitudes proportional to the RCS and are assumed to be the same for 
all the receivers. The latter assumption is consistent with a small network in which the distances between network 
nodes are much smaller than the distances between the nodes and the targets, i.e., d k {t) ^> r l l T . Thus, since they 
are closely spaced, all receive nodes see the same aspect of the target. 
Due to reflection by the target, the Z-th antenna element receives 

zi{t) = f> fc (t-M))+ Cj (*) 

= f> f Xi (t - + *® y +ei{t) , l = Mr (5) 

fc=i i=i 

where ei(t) represents noise, which is assumed to be independent and identically distributed (i.i.d.) Gaussian with 
zero mean and variance a 2 . 

The narrowband assumption on the transmit waveforms allows us to ignore the delay in xi(t), and consider the 

delay in the phase term only. Thus, the received baseband signal at the Z-th antenna can be approximated as 

k M t 

zi{t) w Y, 0k E x^e^^e^- 2 ^^^ 6 ^^ + c,(t) 

k=l i=l 
K 

= Y^e-^^e^^e^^it^ie^ + e^t) (6) 
fe=i 
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where A is the transmitted signal wavelength, fk = 2vkf /c is the Doppler shift caused by the fc-th target, and 

v(fl fc ) = [eW^,...,^^^ (7) 

andx(f) = Mt),..,x Mt W] T (8) 

On letting L denote the number of snapshots and T s the sampling period, the received samples collected during 
the m-th pulse are given by 

zi((m - 1)T + 0T S ) 

' =^ 7te JTtf(^) e ^( m -i)r D(/fc)Xv( ^ ) + eim (9) 



z/,. 



zi{{m-l)T+{L-l)T s ) 



fc=i 



where 



7fc = ^ e -^^(o) ) 

D(/ fe ) = dxag{[e^ 0T ',...,e^M L -^]}, 

e im = [ei((m-l)T + QT s ),...,ei((m-l)T+(L-l)T s )} T , 

andX = [x(0T s ),...,x((i- l)T s )f (L x M t ). 



(10) 



In this paper we assume that the Doppler shift is small, i.e., fkT s << 1 for k = 1, due to slowly moving 

targets. 

III. Compressive Sensing for MIMO Radar 
Let us discretize the angle-Doppler plane on a fine grid: 

a = [(01,61), . . . , (a N ,b N )]. 

We can rewrite (9) as 

N 

zi m = J2 s ne J ^ M e^ K{m - 1)T n(b n )Xv(a n )+e lm 

n=l 

where 

{7fc, if the fc-th target is at (a n , b n ) 
■ 
0, otherwise 

In matrix form we have 

Zfm = *i m S + ei m 

where s — [si, . . . , sn] t and 

*(m = [e j2 f"^ ai ) e j27rbl(m - 1)T D(fe 1 )Xv(a 1 ), . . . , e^ I '' r(oN) e j2,r '' N(m - 1)T D(6 A r)Xv(aAr)]. 



(11) 



(12) 



(13) 



(14) 



(15) 



Assuming that there are only a small number of targets, the positions of targets are sparse in the angle-Doppler 
plane, i.e., s is a sparse vector. Let us measure linear projections of z/ m as 



lm •■ 



(16) 
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where &i m is an M x L (M < L) zero-mean Gaussian random matrix that has small correlation with ^i m , and 
&l m = &i m eim- M must be larger than the number of targets. 

If the Z-th node in the network knows who the transmit nodes are and also knows the transmitters' coordinates 
relative to a fixed point in the network, then the node can construct the matrix *f?i m (15) and recover s via l\- 
optimization based on the node's own received data r; TO (see (16)). Information on other nodes' locations could 
be provided by higher network layers. If no such location information is available to the node, or the interference 
is strong, then the receive node will pass the linear projections r; m to a fusion center, which has global and local 
information. Combining the output of N p pulses at N r receive antennas the fusion center can formulate the equation 

r = [rJ 1 ,...,r^ Np ,...,rJ frl ,...,r% rNp ] T = 0s + E (17) 

where = [(*u*u) T , . .., (*i (JVp -i)*i(iv p -i)) T , ■ ■ ■ , {^ Nr i^N T i) T , ■ ■ ■ , (*jv r (jv p -i)*w r (w p -i)) T ] T andE = 
[e[ 1: . . . , efjy , . . . , ejy- rl , . . . , &N r N ] T - Thus, the fusion center can recover s by applying the Dantzig selector to 
the convex problem of (17) as ([24]) 

s = mm||s||i s.t. \\& H (r - 0s) < \i. (18) 

According to [24], the sparse vector s can be recovered with very high probability if fj, = (l+t~ 1 )y/2 log Na 2 cr max , 
where t is a positive scalar, a max is the maximum norm of columns in the sensing matrix and a 2 is the variance 
of the noise in (17). If &<f> H = I then a 2 = a 2 . Determining the best value of \i requires some experimentation. 
A method that requires an exhaustive search was described in [24]. A lower bound is readily available, i.e., fj, > 
\J 2 log No 2 a m ax ■ Also, fx should not be too large because in that case the trivial solution s = is obtained. Thus, 
we may set \i < ||0 r||oo. 

A. Resolution 

The uniform uncertainty principle (UUP) [11] [12] indicates that if every set of columns with cardinality less than 
the sparsity of the signal of interest of the sensing matrix (0 defined in (17)) are approximately orthogonal, then 
the sparse signal can be exactly recovered with high probability. For a fixed M the correlation of columns of the 
sensing matrix can be reduced if the number of pulses N p and/or the number of receive nodes N r is increased. 
Intuitively, the increase in N p and N r increases the dimension of the sensing matrix columns, thereby rendering 
the columns less similar to each other. A more formal proof is provided in Appendix I. Moreover, increasing the 
number of transmit nodes, i.e., M t , also reduces the correlation of columns; this is also shown in Appendix I. 

In general, to achieve high resolution a fine grid is required. However, for fixed N p , N r and M t , decreasing 
the distance between the grid points would result in more correlated columns in the sensing matrix. Based on the 
above discussion, the column correlation can be reduced by increasing N p , N r or Mt. Also, based on the theory 
of CS, the effects of a higher column correlation can be mitigated by using a larger number of measurements, i.e., 
by increasing M. In particular, it was shown in [11] that M should satisfy M > — Sg^O ; where e denotes the 
maximum mutual coherence between the two columns of the sensing matrix and C is a positive constant. 
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One might tend to think that in order to achieve good resolution one has to involve a lot of measurements, or 
trasnmit/receive antennas, or pulses, which in turn would involve high complexity. However, extensive simulations 
suggest that this is not the case. In fact, the proposed approach can match the resolution that can be achieved with 
conventional methods, while using far fewer received samples, than those used by the conventional methods. 

B. Maximum grid size for the angle-Doppler space 

The grid in the angle-Doppler space must be selected so that the targets that do not fall on the chosen grid points 
can still be captured by the closest grid points. This requires sufficiently high correlation of the signal reflected by 
each target with the columns of corresponding to grid points close to the targets in the angle-Doppler plane. 
However, this requirement goes against the UUP, which requires that every set of columns with cardinality less 
than the sparsity of the signal of interest be approximately orthogonal. Thus, there is a tradeoff of the correlation 
of columns of the sensing matrix and the grid size. 

Absent prior information about the targets, we can determine the maximum spacing of adjacent grids in the angle- 
Doppler space by considering the worst case. Assume that we discretize the angle-Doppler space uniformly with 
the spacing (Aa, Ab) as a = [(01, 61), . . . , (ojv, 6jv)]. The worst case scenario is that the targets fall in the middle 
between two adjacent grid points. Therefore, a practical approach of selecting the grid points is to calculate the 
correlation of columns corresponding to (af + 4f > with the columns corresponding to (04, bi), i = 1, . . . , N. 

This can be done by computing the correlation at lag zero of columns corresponding to (a, + 4p ? h + 4r) with 
the columns corresponding to (a*, bi), for i = 1, . . . , N, and then taking the average. Then, we can vary the step 
(Aa, Ab) until the average correlation reaches some threshold. This threshold should be high enough to capture 
the targets that do not fall on the grid in the angle-Doppler space, and at the same time, it should satisfy the UUP. 
The adoption of such grid points would ensure that the angle-Doppler estimates of targets would always fall on the 
grid of the constructed basis matrix. 

When the targets are between grid points, the increase in N p or N r will not necessarily improve performance. 
However, simulations show that we can obtain very good performances with very small N p and N r . To achieve a 
similar performance, the conventional matched filter method will require much greater N p and N r . 

C. Range of unambiguous speed 

Let us assume that the Doppler shift change over the duration (T) of the pulse is negligible as compared to the 
change between pulses. This is a reasonable assumption given that we have assumed fkT s << 1, k = 1, . . . ,K. 
Given two grid points (a,,6j) and (ai,bj) in the angle-Doppler space, where bi 7^ bj, the corresponding columns 
of \f are different if eJ 2 ' KhiT ^ e j 2vb 3 T _ Let Vi be the speed corresponding to the Doppler frequency bi and 
A^ = Vj - Vi. It holds that 




2A«/T 



^n,n = ±1,±2, . . . 



(19) 



c 
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Therefore, the range of the unambiguous relative speed between two targets that appear at the same angle satisfies 

2A>i fT r 

» 1 < 1 A» < — . (20) 

The selection of T affects the range of the unambiguous speed; the smaller the T the larger the range of the 
unambiguous speed is. We also need a relatively small T to satisfy the assumption that the Doppler shift does 
not change within the duration of the pulse. On the other hand, a larger T is needed to satisfy the narrowband 
assumption about the transmitted waveforms. Therefore, T needs to be chosen to balance the above requirements. 

D. Complexity 

The proposed approach requires solving the convex programming problem of (18). The more targets one would 
hope to be able to detect the higher the complexity would be. Further, the signals involved are complex. In this case 
(18) can be recast as a second-order cone program (SOCP) [13], which requires polynomial time in the dimension 
of the unknown vector. 

The requirement of a fine grid further increases the computational complexity. This problem can be mitigated by 
first performing an initial angle-Doppler estimation using a coarse grid, and then refining the grid points around 
the initial estimate. Restricting the candidate angle-Doppler space reduces the samples in the angle-Doppler space 
that are required for constructing the basis matrix, thus reducing the complexity of the i\ -optimization step. 

In addition to the computation complexity, the receiver for obtaining the required samples is also more complex. 
The schematic diagram of the receiver is shown in Fig. 1 (see also [18]). 

IV. Performance Analysis in the presence of a jammer signal 

In [24], Candes and Tao showed that if the basis matrix obeys the UUP and the signal of interest s is sufficiently 
sparse, then the square estimation error of the Dantzig selector satisfies with very high probability 

N 

II s-s ||| 2 < C 2 2logNx (ct 2 + ^ min(s 2 (i), cr 2 )) (21) 

i 

where C is a constant, N denotes the length of s and a 2 is the variance of the noise. It can be easily seen from (21) 
that an increase in the interference power degrades the performance of the Dantzig selector. Thus, in the presence 
of a jammer that transmits a waveform uncorrelated with the radar transmit waveforms, the performance of the 
proposed CS method will deteriorate. Next, we provide analytical expressions for the signal-to-jammer ratio at the 
receive nodes, and propose a modified measurement matrix to suppress the jammer. 
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A. Analysis of Signal-to-Jammer Ratio 

Suppose that each transmitter transmits N p pulses. In the presence of a jammer at location (d,8) the signal 
received at the Z-th receive antenna can be expressed as 

$ lie j2*/ fc 0T 



rz = 



r/i 



K 

E 

fc=i 



Ike 



$ e j2Trf k (N p ~l)T 



D(/ fc )Xv(fl fc ) 



+ e -i^(<J-»rW) / 3 





+ 











(22) 



where x m = [x m (0T s ), . . . ,x m ((L — 1)T S )] T contains the samples of the signal transmitted by the jammer during 
the m-th pulse, and j3 denotes the square root of the power of the jammer over the duration of one signal pulse. 

We assume that for all m, E{x^ n (i)x m (j)} = 1/L for i = j, and otherwise. Thus, _E{x^x m } = 1. Also, we 
assume that x m , m = 0, . . . , N p are uncorrelated with the main period of the transmitted waveforms. Thus, the 
effect of the jammer signal is similar to that of additive noise. In the following analysis we assume that the jammer 
contribution is much stronger than that of additive noise, and therefore we ignore the third term r/„ on the right 
hand side of (22). Later, in our simulations we will consider additive noise in addition to a jammer signal. 

We assume that all receive nodes use the same random measurement matrix over N p pulses, i.e., €>; = 3>n = 



. . . = * /A r p . Let Af' fc = X H D H (/ fe )*f *jD(/ fc /)X and q- f denote the (ij)-th element of Af fe . Thus, 
the average power of the desirable signal conditioned on the transmitted waveform can be represented by 

K N p -1 

P.(l) = E{rf s r ls \X} = E{ V 7*7^"^ ( V] e - j2 <^- f ^ mT )v H (6 k )Af v(6 k ,)} 



k,k' = l 



pi{k,k>) 



m=0 



Qkk' 



Mfcfc' 



K 



= N p E{Y^ \Pk\ 2 Qkk} +E{J2 Plfc k')nkk>Qkk>} 
fc=i fe^fc' 

where pi(k, k') and Qkk' can be further written as 



pi(k,k') = J*MMo)-d h ,m-(.rSV*)-vWk>))ftp kl 



(23) 

(24) 
(25) 



t It t It \ 

As defined in Section II, the position of the ith transmit or receive (TX/RX) node is denoted by (ry , a/ ) in 
polar coordinates. Thus it holds that 



.fc'fc 

l 3i 



2r/ sm( 2 ) sin(a, 2 > 1 = J 

r/ r cos(6 k < - aij) - r- /r cos(6» fe - a,-) j 



(26) 
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Let -00 be deterministic. Based on the assumed statistics of r-s and a, (see (2)), the distribution of h 



t/r .i 

^— sin(a/ r - ip ) is given by ([26]) 



f h (h) = -^Jl-K 2 , -Kh<\ (27) 

7T 



and 



E (M =2 W_ (28) 
L ' a 

where Ji(-) is the first-order Bessel function of the first kind. Thus, based on (28) we can obtain 

1 i = j and k = k' 

{e j2 ? a -' fc } =e\ tP^^ \ = <j ^(4sin(2i^)) i = j and k ^ k' (29) 

? 2 (2) i^j 



where = 2^p. 

Therefore, the average power of the desirable signal P s (l) taken over the positions of TX/RX nodes can be found 
to be 

P s (l) = N p E ij^l^Qkk] +E J E MMWQfcfc' 

K 

= N p J2\Pk\ 2 E{Q kk }+J2 E{ P i(k,k')}nkk'E{Q kk ,} 

k=l k=£k' 

= ^El^| 2 E^{^- a - }+ E ^^ e ^ (dfc(0) - dfc ' (0)) i?{e^ Q " fc }^E^f E{e^ a " k } 

k=l i,j k^k' i,j 

= ^Ei/3 fc i 2 E<?"+E<^ 2 ( 2 )] 

k— 1 2 

k^k' i i^j 

(30) 

where ftfc' = ^(4sin( e ^'~ efc )). 

For many practical radar systems with wavelength A less than 0.1m, (e.g., most military multimode airborne 
radars), 2-kt/X is a large number if r > 5m. Since the function <;(x) decreases rapidly as x increases, the terms 
multiplied by <; 2 (2) are small enough to be neglected in the above equation. Therefore, (30) can be approximated 
by 

P.(t) « ATpf>| 2 E^+E»' e ^ (3« 

fe=l i fe^fc' i 
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Similarly, the average power of the jammer signal over TR/TX locations is given by 

Pj (l) = Eivfa} = { e-^ d -^0)(e-W-^0T ^ *£*f **» 

m—l 

= |/3| 2 g 5c^$f *,x m . (32) 

m— 1 

The SJR given the node locations is the ratio of the power of the signal to the power of the jammer. Since the 
denominator does not depend on node locations, the average SJR equals SJR= P s (l)/Pj(l). 

Some insight into the above obtained expression will be given in the following for some special cases. 

B. SJR based on a modified measurement matrix 

Since the jammer signal is uncorrelated with the transmitted signal, the SJR can be improved by correlating the 
jammer signal with the transmitted signal. Therefore, we propose a measurement matrix of the form 

= *JX H (M x L) (33) 

where <fr| is an M x M t Gaussian random matrix. Note that <!>/ is also Gaussian. As stated in [12], a random 
measurement matrix with i.i.d. entries, e.g., Gaussian or ±1 random variables, is nearly incoherent with any fixed 
basis matrix. Therefore, the proposed measurement matrix exhibits low coherence with "J//, thus guaranteeing a 
stable solution to (18). Based on (33), the average power of the desirable signal P s (l) is given by (30), except that 
Qkk< is based on A, ' = X ff D ff (/fc)X($J) ff ^ > JX H D(/fc')X. The average power of the jammer signal is given 
by (32) where is replaced by <!?/. 

Let us assume that the Mr transmit nodes emit periodic pulses containing independent quadrature phase shift 
keying (QPSK) symbols, and that X H X = l Mt . Also, we assume that = = l M - 

Let Xi(n) be expressed as §i n /"jL, where $i„ is a random variable with mean zero and variance one. Then the 
average power of the jammer signal Pj(l) can be rewritten as follows: 

N p 
m— 1 

N p L-1 N p L-1 

= I/ 3 ! 2 E E S m(*) £ ™ + |/3| 2 Y E 5 m(*) 5 ™0') C y 

m—l i—j—0 m—l 

Np L-1 Np L-1 

= j)P\ 2 E E ^miCii + Y E VmjCij (34) 

m—l i— m—l 

where Cjj is the (i, j)-th entry of $>f$>i. Since the entries of are i.i.d Gaussian variables with zero means and 
variances j-, ca,i = 1,...,L are i.i.d chi-square random variables with means ^ and variances Cy, i ^ j 
are of mean zero and variance M/L 2 . Let us express Cij,i ^ j as Qijy/M / L, where Qij has zero mean and unit 
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variance. It holds that 



N p r— N p L-l 

Pj(l) = \P\ 2 J2 E W*rni»miCii} + E E Ku^lQii 

m— 1 m— 1 



L L \ L(L - 1) 



i&=0 
m— 1 

iV p |/3| 2 ^ (35) 



where we have used the fact that for large L, 

L-l 



fE^^ £R«<W«} = ^ (36) 

(37) 



and - s E] ^mi^mjQij ~* ^C^ftj} = • 

Using the measurement matrix <I>j in (33) will not affect the average Pj(l) over the jammer signal due to the 
fact that ]n jCil = Tr{X(*;) ff *;X ff } = Tr{X H X(*;) ff *J} = Tr{$;($;) H } = Tr{Ij\/} = M. 

In the following, we will look into the SJR improvement using 3>/ as opposed to for two different cases, 
i.e., stationary targets and moving targets. 

1) Stationary Targets: First, let us consider the SJR using the random measurement matrix 

When the targets are stationary, the Doppler shift is zero and so Af' fe = A, = X H *f *;X. Therefore, the 
average power of the desired signal can be approximated as 

Ps (i) « ^ p f:i^i 2 E^ + ^E^^ eJ * (dfc(0Kdfc ' (0)) ^E^ 08) 

k— 1 i k^k' i 

where qi,j is the (i,j)-th entry of A;. 

Letting denote the z-th column of X, ^ - can be expressed as 

M t M t L 

= Tr { A J = E x ^*^*' x * = E E j; i( m )«('') 

i i—1 i—1 m,n—l 

Mt L M t L 

= EE x ii m ) x i( m ) c m m + E E x i ( m ) x i( n ) c mn- (39) 
i—1 rn— 1 i—1 m^n 

The entries of X have zero means and mutually independent; therefore, for sufficiently long L and M t it holds 
that 

> J qi,i = — 2^ c »»» = 7 ■ ( 4 °) 

?' m— 1 

Based on (40), a concise form of P s (0 is given by 

P<(0 w + ( 4i) 
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where <p = Ek,v,k*k> ^/W^^'^V 

Thus, the SJR corresponding to the random measurement matrix <J>; is 

5Ji?; - w) w — ■ m 

When using the measurement matrix 4>; = <I>jX H , the quantity corresponding to A, ' is 

Af< fc ' = A; = x ff x(*o H $;x ff x = . (43) 

It holds that J^V q,^ = = Tr{*J(*J)' ff } = M. Similarly, the average power of the desired signal 

can be approximated as 

K 

PS) « ^V p Af(^|/3 fc | 2 + ^). (44) 
fc=i 

Therefore, the SJR corresponding to the random measurement matrix is 

SJRi = ps) w — ■ ( } 

From (42) and (45), it can be seen that the use of <E>; instead of can improve SJR by a factor of L/Mt when 
L ^> Mt- The SJR can be improved by an increase in L. However, increasing L will require a higher sampling rate 
when the pulse duration is fixed. It is interesting to note that the SJR of (45) does not depend on the the number 
of measurements, M. 

2) Slowly Moving Targets: For simplicity, we consider only the scenarios in which feT << 1. 
Based on the measurement matrix <J>;, and considering the Doppler shift, we have A k ' k = ~K H T> H (fk)^f 1 &iT>(fk> 
When the normalized Doppler frequency fkT s < 1, we have 



J2<lif = Tr{Af< fc '} = Tr{X ff D ff (/ fc )*F*iD(A,)X} » MMl. (46) 



L 

Thus, P s (l) for the moving targets with f s T << 1 is approximately the same as that of stationary targets 

c fc - 



Let us now consider the measurement matrix 3>/. Let c& denote the (i,j)-th entry of X ff D ff (/fe)X and note 



that c*j is given by = 2~2n=o x i( n ) x j( n ) * e j27r ^ fc ™ Ts . In scenarios in which fkT s << 1 and L is relatively 
large, the following approximations are readily derived: 



1 1— e 



c fe ^ u - _ (4?) 



i l_ e j27r/ fc T s * — J 

»v j 

Since the off-diagonal elements are small compared with the diagonal elements, they can be ignored. 
Then, we obtain the following approximation 

a m = x ff D H (/,)x(*;) H *;x ff D(/ fc )x « (48) 

Therefore, the SJR of moving targets with f s T « 1 is approximately equal to that of stationary targets for both 
random measurement matrices. 
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V. Simulation Results 

The goal of this section is to demonstrate the ability of the proposed MIMO radar approach, denoted in the 
figures as CS, to pick up targets in the presence of noise and/or a jammer, and also show the effect on the various 
parameters involved. In each case the performance is compared against other methods that have been proposed in 
the context of MIMO radar (here referred to as "conventional") in order to quantify weaknesses and advantages. 
For the case of stationary targets, the conventional methods tested here are the methods of Capon, APES, GLRT 
[2] and MUSIC [27], while for moving targets, comparison to the matched filter method [21] is conducted. 

In our simulations we consider a MIMO radar system with the transmit/receive antennas uniformly distributed 
on a disk of radius 10m. The carrier frequency is / = 5GHz and the sampling rate f s = jr = 20MHz. The pulse 
repetition interval is T = 1 /4000s. Each transmit node uses uncorrelated QPSK waveforms. The received signal is 
corrupted by zero mean Gaussian noise. We also consider a jammer that transmits waveforms uncorrelated to the 
signal waveforms. For simulation purposes we take the jamming waveforms to be white Gaussian [28]. The SNR 
is defined as the ratio of power of transmit waveform to that of thermal noise at a receive node. 

A. Stationary Targets 

The presence of a target can be seen in the plot of the magnitude of s obtained by (18). We will refer to this 
vector as target information vector. The location and magnitude of a peak in that plot provides target location and 
RCS magnitude, respectively. The proposed approach results in a clean plot away from the target locations, and 
well distinguished peaks corresponding to the targets. This is a desirable behavior for target detection, as it would 
result in small probability of false alarm. To demonstrate the appearance of the graph we define the peak-to-ripple 
ratio (PRR) metric as follows. For the fc-th target, PRR^ is the ratio of the square amplitude of the DOA estimate 
at the target azimuth angle to the sum of the square amplitude of DOA estimates at other angles except at the 
jammer location, i.e., PRRu = H , — rs — i — f2> where s is the defined in (13), Sfc and Sj denote the elements 

s s — z^i = i\Sk\ ~\ s j\ 

of s corresponding to the location of the fc-th target and the jammer, respectively. A clean plot would yield a high 
PPR, while a plot with a lot of ripples would yield a low PRR. 

A metric that shows the degree to which a jammer is suppressed, namely the peak-to-jammer ratio (PJR), is also 
used here. PJR is defined as the ratio of the average square amplitude of the DOA estimates at the target angles 

— v K - is t i 2 

to the square amplitude of DOA estimates at the jammer, i.e., PJR = K — — ■ Unlike PRR, PJR is averaged 
over all targets. In this way, the jammer is considered to be suppressed only if the peak amplitude at the jammer 
location is much smaller than the peak amplitude at any target location. 

The results that we show represent 1 , 000 Monte Carlo simulations over independent waveforms and noise 
realizations. To better show the statistical behavior of the methods we plot the cumulative density function (CDF) 
of PPR and PJR, i.e., Probability (PPR < x) and Probability (PJR < x), where PPR is the union of PRR k , k = 
1,..., K. 

1) Targets falling on the grid: We consider the following scenario. Two targets are located at angles 0\ =0.2° 
and 9 2 = —0.2°. The corresponding reflection coefficients are (3% = (3% = 1- A jammer is located at angle 7° and 
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transmits an unknown zero-mean Gaussian random waveform with variance 1 = 400. Additive white Gaussian 
noise is added at the receive nodes. The ratio of the power of transmitted waveforms at each transmit node to the 
variance of the additive Gaussian noise is set to dB. The number of transmit antennas is fixed at M t = 30. For 
the purpose of reducing computation time, the angle space is taken to be [—8°, 8°], and is sampled with increments 
of 0.2° from —8° to 8°, i.e., a = [—8°, —7.8°, . . . , 7.8°, 8°]. The received signal in a single pulse is sampled, and 
M = 30 random measurements of one pulse are used to feed the Dantzig selector. Since the MUSIC method 
requires the number of receive antennas to be greater than the number of targets, when only one receive antenna 
is used we compare the proposed CS method with only the Capon, APES and GLRT methods. The comparison 
methods are using L = 512 samples to obtain their estimates, while the proposed approach uses M = 30 samples. 
The result of one realization for the case of one receive node is shown in Fig. 2. One can observe the cleaner 
appearance of the graph corresponding to the proposed approach, where the two targets appear correctly except 
with a small error in the magnitude of the target RCS. The CDF of the corresponding PRR and PJR are also shown 
in the same figure. One can clearly see that with one receive antenna the comparison methods yield PRR close to 
1, which is indicative of severe ripples. 

In general, an increase in the number of transmit snapshots L leads to improved PRR and PJR for all methods. 
In the following results we fix L to 512. For the comparison methods, L represents that number of samples needed 
to obtain target information. For CS, the number of samples used to extract target information is M. 

For the scenario of Fig. 2, the effect of the threshold /j, is evaluated in terms of the empirical CDF of the PRR 
and the amplitude estimate of RCS, and the results are shown in Fig. 3. One can can see that the increase in fj, can 
lead to fewer ripples but at the same time it degrades the amplitude estimate of RCS. In the following, the values 
of fi used in each case will be shown on the figures. 

For the same target and jammer configuration as above, we now examine the effect of different levels of jammer 
strength. We consider the scenario where iV r = 10 receive nodes participate in the estimation. For the case of CS, 
each node sends to the fusion center M = 30 received samples, while for the comparison methods, each node 
sends to the fusion center L = 512 received samples. In Fig. 4 we show the CDF of PPR and PJR corresponding 
to jammer variance (3 2 — 400, 1600 and 3600 and SNR equal to dB. One can see that for CS, the probability of 
low PRR and PJR increases when the jammer becomes stronger. In particular, there is some non-zero probability 
that the PRR will be close to 10 -7 . Such cases are rare and occur when one of the two targets is missed. The 
increase in the threshold /i can improve the DOA estimates at the target locations and reduce the probability of 
missing one of the targets. The cost, however, would be an increase in ripples. The performance of the proposed 
approach can be improved, i.e., the rare low PPR values can be completely avoided by increasing N r , or M, This is 
demonstrated in Fig. 5, where the strong jammer case of Fig. 4 is considered, i.e., (3 2 = 3600, and N r is increased 
to 30. We should note here that it does not help to increase M beyond M t as the maximal rank of $J is M t . 

Next, we consider the same scenario as above but let the two targets be at variable distance d in the angle 
domain. Figure 6 demonstrates performances for the cases d = 0.2°, 0.3°, 0.4° in the presence of a strong jammer 
with variance [3 2 = 3600. The SNR is dB, iV r = 10 and M = 30. One can see that the comparison methods 
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produce good level PRR. Regarding the PJR, as expected, MUSIC fails, Capon and APES results is PRR« 1 most 
of the time, while GLRT performs well all the time. The proposed CS approach performs well with a few exceptions 
in which a PRR or PJR less than 1 is obtained with very small probability. Again, the CS method performance can 
be improved by increasing N r and/or M. 

Based on the above results, the performance of the proposed approach for the jammer dominated scenario can be 
made at least comparable to that of the conventional methods while using about 5.8% (= 30/512) of the number 
of samples required by the conventional methods. 

Next, we study a thermal noise dominated case, i.e., SNR=— 40dB. Figure 7 shows PRR and PJR performance for 
different values of jammer variance, i.e., 1 = 400, 1600 and 3600. In all cases it was taken N r = 10, M t = M = 30 
and the targets were separated by d = 0.4°. CS yields good performance even in the presence of both a strong 
jammer and thermal noise. The PRR performance of other methods appear to deteriorate at this noise level. The 
performances for targets with spacing d = 0.2°, 0.3° and 0.4° are given in Fig. 8 for N r = 20, M t = M = 30 
and ft 2 = 400. Like in the case of a strong jammer, the decrease in the spacing d does not affect the performance 
significantly. In this thermal noise dominated case, CS appears to perform very well in terms of PRR, and PJR, 
while the comparison methods appear to be very noisy. To further examine this case, we consider two additional 
performance measures, i.e., mean squared error (MSE) and probability of false alarm (PFA), which are computed 
based on the obtained estimate s as follows. A new vector, s* is formed; if is greater than some threshold then 
s\ = 1, otherwise, s\ = 0. The MSE is calculated as MSE = ||s t — S( HI/AT, where s t is an N x 1 vector which 
contains zeros everywhere except at angles corresponding to target locations, where it is 1. The PFA measures 
the probability of 1 occurring in s at non-target locations. Figure 9 shows the MSE based on 8, 000 Monte Carlo 
simulations. Note that the performance of MUSIC is not shown here since MUSIC always yields a peak at the 
jammer location. One can see that the simple thresholding described above helps the comparison methods, and 
if the threshold is picked appropriately all methods can produce a low angle MSE and PFA. However, the MSE 
corresponding to the CS method is less sensitive to the particular threshold than other methods. For the milder 
jammer case (/? = 20), the CS approach exhibits slightly better "best MSE performance" than the comparison 
methods, while in the stronger jammer case ((3 — 60) GLRT outperforms CS for most thresholds. For the strong 
jammer case, the MSE and PFA of CS are compared to those of GLRT for different number of samples, L in 
Fig. 10. One can see that for the strong jammer case ((3 = 60) CS performs comparably to GLRT with L = 256. 
Thus, in the strong jammer case, CS still achieves good performance with fewer samples than GLRT, except that 
the savings in terms of number of samples is smaller. For CS, the trend of an increasing MSE as the threshold 
increases can be explained by the fact that one of the two targets can be missed as the threshold increases. GLRT 
relies on the Gaussian assumption for the noise and jammer signals, which is totally valid in our simulations. Thus, 
unlike the other methods, GLRT can suppress the jammer completely. We should note that the specific values of 
MSE and PFA depend on the kind of thresholding performed. For example, applying thresholding on a nonlinear 
transformation of the estimated vector can give different MSE and PFA, and the best results for each method are 
not necessarily obtained based on the same non-linear transformation. Determining the best thresholding method is 
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outside the scope of this paper. 

2) Targets falling off the grid points: In this section, we consider scenarios in which targets do not fall on 
the grid points. This is a case of practical interest, as the target locations are unknown, thus the best grid in not 
known in advance. We first select the proper step to discretize the angle space following the procedures described in 
Section III-B. The angle space is sampled by increments of 0.2° from -8° to 8°, i.e., a = [-8°, -7.8°, . . . , 7.8°, 8°]. 
Assume that four targets of interest are located at 6k = { — 1.1°, —0.3°, 0.3°, 1.1°}. Their reflection coefficients are 
{(3k = 1) k = 1,2,3, 4}. A jammer is still located at 7°. Since the targets are located between the grid points, we 
cannot plot PRR and PJR as in the case of targets onto the grid points. Therefore, we show the mean plus and 
minus one standard deviation (std) for the amplitude of DOA estimate at each grid point. The results are shown in 
Fig. 11. The power of the jammer was 400 (left column of Fig. 11) and 3600 (right column). Based on Fig. 11, it 
can be seen that with the proper grid points, the proposed method can capture well the targets that do not fall on 
grid points. The next best method is the GLRT which captures the targets but exhibits high variance as indicated 
by the shaded region around the mean. 

B. Moving Targets 

We continue to consider orthogonal QPSK waveforms and a jammer located at 7° with the power 400. The SNR 
is still set to be dB and each receive node collects M = 30 measurements. Figures 12 and 13 show the target 
scene of the proposed CS method and the matched filter approach [21] for targets on the grid points and off the grid 
points, respectively. The matched filter correlates the receive signal with the transmit signal distorted by different 
Doppler shifts and steering vectors. 

1) Targets falling onto the grid points: We assume the presence of three targets located at {6k = —1°, 0°, 1°} 
that are moving at the speed of {vk = 60m/ s, 70m/ s, 80m/s}, respectively. We sample the angle-Doppler space 
by the increment (0.5°, 5m) as 

a = [(-8°, 50m/ s), (-7.5°, 50m/ s), (8°, 50m/ s), (-8°, 55m/ s) . . . , (8°, 55m/ s), (8°, U0m/s)} (49) 

Figure 12 shows the target scene for one realization corresponding to Ni = 1 receive nodes (left column of the 
figure), and also N r = 10 (right column of the figure). We can see that the performance of the match filtering 
method is inferior to that of the CS approach even when using the data of 30 pulses. The proposed CS approach 
can yield the desired performances even with a single receive node and as low as 5 pulses. Comparing the left 
column and right column of Fig. 12, one can see the effect of the number of receive antennas N r . The increase in 
N r can reduce the number of pulses required to produce good performance. 

2) Targets falling off the grid points: In this section, we consider the scenarios in which targets that do not fall 
on grid points. From simulations (the corresponding figure is not given here because of space limitations), we found 
that the column correlation is more sensitive to the angle step than the speed step, since fT s << 1. This indicates 
that in the initial estimation, the grid points should be closely spaced in the angle axis and relatively sparser in the 
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speed axis. Then the resolution of target detection can be improved by taking denser samples of the angle-Doppler 
space around the initial angle-Doppler estimate. 

Like the scenarios with the stationary targets, the angle dimension is sampled by increments of 0.2° and the step of 
the speed dimension is set to 5m/ s. Three targets are moving at the speed of {i>k = 62.5m/ s, 72.5m/ s, 82.5m/s} 
in the direction of {9k = —1-1°, 0.1°, 1.1°}. Fig. 13 demonstrates that the proposed method can capture the targets 
which fall out of the grid points in both angle and speed dimensions and it can outperform the conventional matched 
filter method. Moreover, we can see that the increase in N p or N r will not necessarily improve performance for 
the targets between grid points. This is because an increase in the dimension of the basis vectors will decrease the 
correlation of columns in the basis matrix, which contradicts the requirement for capturing the targets out of the 
grid points III-B. The performance in the case of closer spaced targets, i.e., d = 0.4° is shown in Fig. 14. 

VI. Conclusions 

We have proposed a MIMO radar system that can be implemented by a small-sized wireless network. Network 
nodes serve as transmitters or receivers. Transmit nodes transmit uncorrected waveforms. Each receive node applies 
compressive sampling to the received signal to obtain a small number of samples, which the node subsequently 
forwards to a fusion center. Assuming that the targets are sparsely located in the angle-Doppler space, the fusion 
center formulates an l\ -optimization problem, the solution of which yields target angle and Doppler information. 
For the stationary case, the performance of the proposed approach was compared to that of conventional approaches 
that have been proposed in the context of MIMO radar. The comparison scenario assumed that each receive node 
forwards the received signal to a fusion center, where Capon, APES, GLRT or MUSIC is implemented to obtain 
target information. The proposed approach can extract target information based on a small number of measurements 
from one of more receive nodes. In particular, for a mild jammer, the proposed method has been shown to be at 
least as good as the Capon, APES, GLRT and MUSIC techniques while using a significantly smaller number of 
samples. In the case of strong thermal noise and strong jammer, the proposed method performs slightly worse than 
the GLRT method. In that case, its performance is still acceptable, especially if one takes into account the fact that 
it uses significantly fewer samples than GLRT. For the case of moving targets, the proposed approach was compared 
to conventional matched filtering, and was shown to perform better in both single and multiple receive nodes cases. 
An important feature of the proposed approach is energy savings. If the fusion center implemented the proposed 
CS approach, it would require nodes to forward M samples each, as opposed to L samples that would be needed 
if the fusion center implemented the conventional methods. In order to meet a certain performance, M is typically 
significantly smaller than L, i.e., fewer samples would be needed for the CS implementation as compared to the 
implementation of conventional methods. This translates to energy savings during the transmission of the samples 
from the receive nodes to the fusion center. The obtained savings would be significant in prolonging the life of the 
wireless network. Future work includes extension to extracting range information, and also studying scenarios of 
widely separated antennas and wideband radar signals. The proposed approach assumes that nodes are synchronized 
and the fusion center has perfect node location information. The effects of localization and synchronization errors 
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and ways to mitigate them need to be further studied. 
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Appendix I 

THE EFFECTS OF N r , N p , Mt ON THE CORRELATION OF COLUMNS IN THE SENSING MATRIX 
A. The effect of the number of pulses on the column correlation in the sensing matrix 
The sensing matrix for the Z-th receive antenna &i is given by 



* i * i(A r p _ 1) 



(50) 



where \^ TO , m = 0, . . . , N p — 1, is defined in (15). 

On letting gfe denote the i-th column of 0;, the correlation of columns gfe and gfe/ equals 

N p \v H (a k )B? k v(a k )\ k = k' 
s^ ibk - K ,) T) ^ H (a k )Bfy(a k ,)\ k + k< 
where Bf fe ' = X H T> H (b k )<Pf $/D(6 fe /)X. 

For a given pair (k, k'), k 7^ k', the ratio of | < gfe, gfe > | to | < gfe, gfe' > |, i.e., hkk', reveals the effect of N p 



Pkk' 



< gfe, gfe' > 



I sin(ir(b k -b k ,)N p T) 



(51) 



on the correlation of the two columns. It holds that 

hkk' oc 



N„ 



(52) 



sm(Tr(b k -b k ,)N p T)\ ■ 

Let assume that T has been fixed. As long as (b k — b k /)N p T < 1, hkk' increases with N p , and attains the maximum 
value when (b k — b k >)N p T = 1, because the cross correlation of gfe and gfe' becomes zero. Therefore, the increase 
in N p can improve the performance of CS estimation of (18) as long as (b k — b k >)N p T < 1. This indicates that if 
(bk — b k ')N p T < 1 for each pair of (k, k'),k 7^ k', the increase in N p can always improve the performances of CS 
estimation. For a conventional radar, the number of pulses can also improve the resolution of Doppler estimates 
since the Doppler shift creates greater change between pulses. 
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B. The effect of the number of receive antennas on the column correlation in the sensing matrix 

Next, we investigate the effect of the number of receive antennas N r on the correlation of columns in the sensing 
matrix. For simplicity, we assume only the received data collected during the n-th pulse is considered and the 
random measurement matrix $ is constant over receive antennas. Then the sensing matrix 6 can be represented as 



e 



(53) 



Thus, the correlation of columns g^ and g 3 equals 



Pij = \<Si,Sj> 



N r 

E- 

2 = 1 



jtt w( a j)-»7r( o 0) 



N r n 



J2'K(n-l)T(b j -b i ) v H 



v H (o,)X ff D iJ (i ! )$ ff $D(6 3 )Xv(o 3 ) 



A, v";</,)B'-'v(« ; 
££r e J * W (<* )"<(«*)) | | v ff (oijB* J 1 v(a.; 
where B lJ = X ff D H (fe i )* H *D(^)X. 

Then the ratio of | < gj, gj > | to | < g 4 , g 4 > | is 

1 



i = 3 
i =h 3 



(54) 



^ oc - 



1=1 



J 2 ? {Vi(°-j)-Vi(<*i)) 



(55) 



Since the receive nodes are randomly and independently distributed, -^-| J2i=i e J "f" (''F ( a j ) — )) | approaches 
as N r becomes large. Therefore, the correlation of two columns in the sensing matrix can be reduced when the 
number of receive antennas is increased. 

C. The effect of the number of transmit antennas on the column correlation in the sensing matrix 

Finally, let us see the effect of the number of transmit nodes on the correlation of columns. For simplicity, we 
assume N r = N p = 1, Then v i? (aj)B l, - 7 v(a :! -) can be rewritten as 

v"i«.jTV-'vi«,: = ]T v k {a 3 )vl{ ai )B^ p /L + Maj)v*Mx k (q)x* k (p)B^ q 

k,p k p^q 



J2 v k( a j)v* k '(ai)x k (q)x* k ,(p)B^ q 

k=£k' p,q 



(56) 



MM t 
L 



UK 



* = 3 
i +3 



L ^ u l 

where Vk and X Pi q denote the k-t\\ entry of v and the (p, <7)-th entry of X, respectively. 
Thus, the ratio of | < gi, gj > | to | < g 4 , g 4 > | is 



L 



MMt 
L 



AI 12k v k( a j) v k( a i) 
M t L 



M t ^ M t 



M 
L 



M t ^ M t 



(57) 



(58) 
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It can easily be seen that the numerator approaches as M t approaches infinity. Therefore, the correlation of two 
columns of the sensing matrix can be reduced by employing a large number of transmit nodes M t . 
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* 0,(/) is the fth row of <J>, of length L 



Fig. 1. Schematic diagram of the receiver. denotes the measurement matrix for the ith receive node. 




Fig. 2. One realization of the DOA estimates (left column) and CDF of PPR and PJR (right column). 7V r = 1, M t = M = 30, 1 = 400, 
SNR= dB and fi = 26. 
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Fig. 3. CDF of PRR (top) and amplitude estimate of RCS (bottom). 
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Fig. 4. CDF of PRR and SJR for /3 = 20, 40, 60 with N r = 10, M t = M = 30 and SNR= dB. The corresponding thresholds are 
At = 120, 190 and 280. 
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Fig. 5. CDF of PRR and SJR for /3 2 = 3600 and SNR=0 dB. Two cases are shown, (N r = 10, M t = 30, M = 30) and (N r = 30, M t = 
30, M = 30). The corresponding thresholds are /i = 280 and 800. 
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Fig. 7. CDF of PRR and PJR for = 20, 40 and 60 with N r = 20, Mt = M = 30 and SNR=-40 dB. The corresponding thresholds are 
At = 350, 440 and 550. 
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Peak-to-jammer ratiio(PJR) 



Fig. 8. CDF of PRR and PJR for two targets with the spacing d = 0.4°, 0.3° and 0.2°. N r = 20, Mt = M = 30, SNR=-40 dB and 
/3 = 20. 
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(a)N =20,M =M=30,SNR=-40 dB,p=20,n=350 (a)N=20,M =M=30,SNR=-40 dB,p=20,n=350 
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Fig. 9. MSE of target information vector and probability of false alarm (PFA) for two targets with spacing d = 0.4° for N r = 20, Mt = 
M = 30 and SNR= -40 dB. 
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Fig. 10. MSE of target information vector and probability of false alarm (PFA) for two targets with spacing d = 0.4° for N r = 20, Mt = 
M = 30 and SNR= —40 dB. The number of transmit waveforms and receive samples per pulse for CS is 512 and 30, respectively. 
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Fig. 11. Modulus of DOA estimates for four targets that do not fall on grid points. The dotted line is the mean of DOA estimates. The yellow 
region is the area bounded by the curves mean ± std. 
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Fig. 12. Angle-Doppler estimates for three targets on the grid points. The three targets are located at {-1°, 0°, 1°}. Mt = M = 30, SNR= 
dB and /3 2 = 400. 
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Fig. 13. Angle-Doppler estimates for three targets that do not fall on the grid points. The three targets are located at {-1.1°, 0.1°, 1.1°}. 
M t = M = 30, f3 2 = 400 and SNR= dB. 
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Fig. 14. Angle-Doppler estimates for three targets on and off grid points. N r = 10, M t = M = 30, SNR= dB, /3 2 = 400 and d = 0.4°. 
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